LOW FREQUENCY FLICKERING OF TT ARIETIS:HARD AND SOFT 
X-RAY EMISSION REGION 



ALTAN BAYKAL and UMIT KIZILOGLU 

Physics Department, Middle East Technical University, Ankara 06531, Turkey 

Abstract. Using archival ASCA observations of TT Arietis, X-ray energy spectra and power spectra 
... of tiie intensity time series are presented for tfie first time. The energy spectra are well-fitted by a two 

' continuum plasma emission model with temperatures ~ 1 keV and ~ 10 keV. A coherent feature at 

0^ ' ~ 0.643 mHz appeared in the power spectra during the observation. 
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O ' 1- Introduction 

TT Arietis has been categorized as a nova-like variable from photometric obser- 

^ . vations (Smak & Stepien 1969, Cowley et al., 1975). The system's spectroscop- 

in ! ic orbital period is found to be 0.13755114(13)'' (Thorstensen et al, 1985) from 

^ ■ measurements of the radial velocities of emission lines. This is longer than the 

^ . photometric period of 0.1329'^ (Tremko 1992). The beat frequency of the photo- 

O I metric and spectroscopic orbital periods is ~ 0.25 d-i which is very similar to the 

■ intermediate polar TV Col (Hellier et al. 1991). This led Jameson et al. (1982) to 
^ I the suggestion that TT Ari is an intermediate polar. The optical brightness of TT 

■ Ari varies irregularly on long time scales. It can remain at visual magnitudes ~ 10 
(high state) for many years, interrupted irregularly by low states with visual magni- 

O I tudes ~ 16 which last less than a year (Shatter et al. 1985). In the high state, TT Ari 

• has shown flickering activity (or quasi-periodic oscillations) with periods between 

! 14 min and 27 min (Semeniuk et al., 1987; Hollander & van Paradijs 1992). 

TT Ari was found to be a hard X-ray source by the Einstein satellite (Cordova 
et al., 1981). Jensen et al. (1983) investigated the correlation of X-ray and optical 
flickering activity using simultaneous observations with the Einstein X-ray obser- 

■ vatory and the Mount Wilson Optical Telescope. They found that the flickering 
activity of X-rays at ~ 17 min (~ 1 mHz) is delayed by ~ 1 min with respect to 
the optical flickering activity. They proposed that the optical flickering in TT Ari 
is produced in the inner accretion disk, and a fraction of the energy is transported 
to a different region where the X-ray flickering is produced. The 1 min time delay 
between the optical and X-ray flickering would be the time required to transport 
energy from the optical flickering region to the X-ray emitting region (Jensen et 
al., 1983). 

In the present work, we analyze archival ASCA data to construct the X-ray 
egergy spectra and power spectra of TT Ari. Section 2 describes the observations. 
Section 3 describes the X-ray spectra and X-ray timing. 



2 



ALTAN BAYKAL AND UMIT KIZILOGLU 



xronos.q Ic 1 



9372.9 9373 9373.1 

Time (days) 



9373.2 9373.3 9373.4 



Figure 1. The light curve of TT Ari with SISO count rates binned in intervals of 86 sec. The times 
indicates MJD 40000+. 



2. Observations 



TT Ari was observed with ASCA on January 20 to 21, 1994 with an effective 
exposure time 14 ksec. The ASCA instrumentation (Tanaka et al. 1994) con- 
sists of four imaging telescopes, each with a dedicated spectrometer. There are two 
soUd-state imaging spectrometers (SIS), each consisting of 4 CCD chips, giving an 
energy resolution of 60-120 eV across the 0.4-10 keV band. The SIS with 2-CCD 
mode and GIS data were taken with a time resolution 8 sec and 62.5 msec respec- 
tively. Two gas scintillation proportional counter imaging spectrometers, GIS, have 
an energy resolution of 200-600 eV over the 0.8-10 keV band. Data were extracted 
within a ~ 4 arcmin radius region for each SIS and within a ~ 6 arcmin region 
for each GIS. The typical mean count rates of TT Ari, for the SIS and GIS are 
~ 0.59 and ~ 0.48 count sec~^, respectively. Standard cleaning for ASCA data 
was applied to eliminate X-ray contamination from the bright Earth, effects due 
to high particle background, and hot flickering SIS pixels. The reduction of the 
ASCA archival data and the correction of the photon arrival times to the Solar Sys- 
tem bary center was performed using XSPEC, XRONOS, XIMAGE and XSELECT 
softwares. In Fig. 1, we present the light curve of the observation. 
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3. Results 

3.1. X-Ray Spectra 

The large number of counts obtained from ASCA observations and the broad ener- 
gy range (0.5-10 keV) allows for fits using various spectral models. The spectrum 
was fitted with a two component Raymond Smith (1977) emission model as expect- 
ed from radiatively cooling shock regions. The measured unabsorbed flux during 
the observation was ~ 9.74 x 10~^^ erg sec~^ cm~^ in the band of 0.5-10 keV. 

In order to see the lower energy part of the spectrum (< 2 keV) more clearly and 
to estimate the column density Nh more accurately, ROSAT archival observations 
of TT Ari were extracted and a two component Raymond-Smith model was fitted 
to combined ASCA and ROSAT data. Our preUminary analysis showed that both 
data sets have similar spectra. In these fits, we assume that the spectral parameters 
of ASCA and ROSAT observations were the same up to a normalization factor. 
Results of the fit, shown in Table 1, imply that the ASCA flux was 1.1 ± 0.02 
of the ROSAT flux. The column density of ~ (4.13 ± 0.11) x lO^o cm'^ 
obtained from the joint fit of the ROSAT and ASCA data sets is consistent with the 
previously deduced value from the ROSAT observation alone (Baykal et al., 1995). 
The two component Raymond Smith model fits the data well with xt ~ 1-03, 
giving the plasma temperatures of about ~ 1 keV and ~ 10 keV, for the two 
components. 

Keeping in mind that the spectral characteristics may change in time and while 
the column density probably remains constant; the column density deduced from 
the joint ROSAT-PSPC and ASCA-S1S0+G1S2 data was used for further spectral 
work with the ASCA observations. In Fig. 2, we present the energy spectra fits with 
the two component model. Only SIS0+GIS2 data are used for clarity. 

Table 1 gives the spectral parameters for the best fit to the S1S0-I-G1S2 data to 
compare with the joint fit and two component Raymond Smith model to SISO+SIS 1+G1S2-I-G1S3 
detectors. 

Spectral results from a two temperature Mewe-Kaastra model are also given. 
This model gives slightly lower plasma temperatures, but it does not improve the 
goodness of fit. 

The residuals do not indicate an extra fine feature at 6.67±0.03 keV. Inset in Fig. 
2 shows data and the model fitted to it. To check whether a line feature component 
is evident for 6.43 keV which may indicate the presence of of emission from the 
white dwarf surface, a variable Raymond Smith model with zero Fe abundance and 
a Gaussian Une feature is tested. The data is again well fitted with a Gaussian fine 
at 6.67 ± 0.03 keV of a one sigma width 68 ± 39 eV and a line flux (2.42 ± 0.10) x 
10~^ photons cm~^sec~^. 

In order to see the spectral modulation over the orbit, hardness ratios (2-10/0.2- 
2 keV) are plotted in Fig. 3. The data are folded at the spectroscopic orbital period 
P=0. 137551 14(13)'^ with an epoch JD= 2,443,729.0175 (Thorstensen et al., 1985). 
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Figure 2. The average ASCA/ SISO, GIS2 spectra of TT Ari. The upper panel shows the fit of the 
two component Raymond-Smith model. Data is binned for 5a confidence. The lower panel is the 
residuals from the fit. Inset shows 6.67 keV feature, binned for 3a confidence, together with the 
model fitted to it. 



Table I 

Spectral Fit Parameters for Two Temperature Models " 



model 


kTi(keV) 


kT2(keV) 


Ai(10-Y 


A2(10"^/ 




) xl 


R-S'' 


1.03 ± 0.05 


10.99 ± 1.17 


4.92 ± 1.36 


6.83 ± 0.14 


4.13 ±0.11 


1.03 


R-S = 


0.99 ±0.12 


9.29 ± 1.02 


3.70 ±2.22 


7.16 ±0.19 


4.13 fixed 


0.96 


R-S 


0.97 ±0.08 


8.49 ± 0.63 


3.67 ± 1.51 


7.26 ±0.13 


4.13 fixed 


1.36 


Meka'^ 


0.76 ± 0.05 


8.10 ±0.49 


3.37 ±0.91 


7.32 ± 0.09 


4.13 fixed 


1.48 



" Spectral fits were performed using the program XSPEC. 

Raymond-Smith model for PSPC-1-SIS1-1-GIS2, SIS1-1-GIS2 and SIS0-i-SISl-i-GIS2-i-GIS3 
detectors respectively. 

" Mewe-Kaastra model for SIS0-i-SISl-i-GIS2-i-GIS3 detectors. 

^ Ai and A2 are the emission measures in units of 10~^* /{AttD"^) J n^dV, where D is the dis- 
tance 

to the source (cm) and n is the electron density (cm~^). 



As seen from the Fig. 3, there is no significant variation of the hardness ratio. The 
binary orbital phases 0.75-0.9 were softer in ROSAT observations (see Baykal et 
al., 1995). However, as these orbital phases are occulted in the ASCA observation 
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Figure 3. The upper panel is the normalized light curve folded at the spectroscopic orbital period 
P=0.13755114(13)'* with the epoch JD= 2443729.0175 (Thorstensen et al., 1985). The lower panel 
is the hardness ratio (2-10 keV/0.2-2 keV) of the light curve. 



due to the satellite orbit around the earth,so we can not address this issue with the 
ASCA data. 

3.2. X-Ray TIMING 

The ASCA observations were interrupted by a number of gaps. In order to exam- 
ine the high frequency oscillations we calculated Fast Fourier transforms for each 
observation window and then averaged them. 

The length of each stretch was chosen as 2048 sec with 4096 bins. The aver- 
age power density spectrum is normalized according to Leahy et al. (1983) and 
rebinned logarithmically in frequency (see Fig. 4). In the power spectrum, the 
counting statistics noise (or white noise) level is around 2 since there are 2 degrees 
of freedom for each power estimate and each estimate obeys a distribution (see 
van der Klis 1989). For frequencies higher than 10 mHz, the variability of TT Ari 
dissolves into the noise level which is consistent with the ROSAT observations 
(Baykal et al. 1995). For frequencies lower than 10 mHz, a red noise component 
(or flicker noise) exists. 

In order to examine the flickering activity and the oscillations at lower fre- 
quencies, discrete power spectra (Deeming 1975) were estimated and normalized 
according to Leahy et al. (1983). 

A power spectrum with 16 sec data bins is presented in Fig. 5 which has a mean 
value of 2 at frequencies higher than 10 mHz. A red noise component (or flicker 
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Figure 4. The average power spectrum of TT Ari. Data were binned 0.5 sec for each GIS2 and 
GIS3 detectors. The power spectra were averaged for each observational window and presented by 
log-log scale (horizontal line represents the counting statistics noise). The windows were of length 
2048 seconds. 
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Figure 5. The average power spectrum of TT Ari for SISO, SISl, GIS2 and GIS3 detectors. Data 
were binned in 16 sec intervals for each detector. Horizontal line represents the counting statistics 
noise. Inset shows the lower frequency part with a Unear scale. 
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Figure 6. The upper panel is the normalized light curve folded at 1553 seconds. The lower panel is 
the hardness ratio (2-10 keV/0.2-2 keV) of the ligh curve. 

noise) dominates at frequencies between ImHz and 10 mHz. The noise is saturated 
(secondary white noise ) at frequencies lower than ImHz (see Fig. 5). 

In order to search for an oscillation , we fitted a constant line to the power at 
frequencies lower than 1 mHz yielding the mean power <Power> = 22. In order 
to see the probability of detection of a significant signal above this level, an expo- 
nential probability distribution, Prob(P)=(l/2o"^) exp (-P/2(T^), which is essentially 
a distribution for two degrees of freedom (Scargle 1982) where <Power>=2(T^ 
=22, is applied. For given parameters, the confidence level of the signal for the 
maximum power (Pjnax=^'^^) at 0.643 mHZ is 1-Np exp (-Pma^/<Power>) ~ 0.9 
(Leahy et al., 1983), where Np=40 is the number of iteration frequencies between 
1 X 10~^ and 1 x 10~^. This confidence level is close to ~ 2a signal detection and 
the error in the period estimate is APeriod = Period^/ Tq^s ~ 170 sec where T is 
the time span of the observations. The symmetric side lobes of the peak frequen- 
cy are the convolved window function of the observation gaps. These are seen at 
ferequencies 0.643 ± 0.17 mHz, where 0.17~^ mHz~^ ~ 5582 sec is close to the 
orbital period of the ASCA satelUte around the earth (~ 5700 sec). We also decon- 
volved the window function by making use of the "CLEAN" algorithm (Roberts et 
al, 1987) and verified the 0.643 mHz oscillation. In Fig. 6, the light curve folded at 
this frequency is presented together with its hardness ratio (2-10 keV/0.2-2 keV). 



4. Discussion 

We find that TT Ari is one of the more "sophisticated" and interesting sources 
among the catalysmic binaries. Our basic results are summarised below. 
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A soft X-ray excess around ~ 1 keV has been seen In magnetic catalysmics 
binaries (Singh & Swank 1993, Ishida et al, 1994, Mukai et al., 1994).These sys- 
tems have spin periods 10 — 10^ sec. According to the theories of magnetic cata- 
clysmic binaries (such as DQ Her systems), most of the X-ray emission is received 
from the shock regions of the accretion columns near the white dwarf magnetic 
poles (Patterson 1994). The cooling plasma in the accretion column has both soft 
~ 1 keV and hard ~ 9 keV X-ray components (Singh k. Swank 1993, Ishida et 
al., 1994). Some spectra show the iron fluorescence Une at 6.43 keV which is an 
indicator of emission from the white dwarf surface (Ishida et al., 1994). Although 
the X-ray spectra of TT Ari shows both soft and hard X-ray components, it deviates 
from the category of magnetic cataclysmic binaries for at least two reasons; i) the 
spin period of the source is not observed, ii) it does not show a 6.43 keV emission 
line. 

EINSTEIN and ROSAT observations have shown that the upper limits of the 
X-ray luminosity from the late type seconderies are ~ 10^^ erg/sec (Fleming et 
al., 1989, Hempelmann et al., 1995). This is much less than the X-ray luminosity 
of TT Ari (L X ~ 1.8 X 10"^^ erg/sec for a distance of 125 pc). This rules out the 
possible explanation that the X-ray flux comes from the corona of the secondary 
star. 

According to Jensen et al, (1983) the X-ray flickering lags the optical flickering 
by ~ 1 minute. This condition implies the possibility of X-ray emission regions 
above and below the accretion disk. Indeed, the 6.67 keV emission line suggests 
the possibility of emission from an accretion disk corona as observed from several 
X-Ray binaries (Kallman & White 1989). 

The observed quasi-periodic oscillations in the optical band show correlation 
with mass accretion rates (Hollander & Paradijs 1993) which can be explained 
by a beat frequency model, as used to explain QPOs in Low Mass X-Ray Bina- 
ries (Alpar &: Shaham 1985). However, this correlation was not seen in the X-ray 
band. The X-ray flux is 10% higher in the ASCA observations compared to the 
ROSAT observations while the peak frequency of flickering oscillations (or quasi- 
periodic oscillation) decrease by 36% (see also Baykal et al., 1995). This is an 
another indicator that the X-ray flux is not directly correlated with the mass accre- 
tion rate. Furthermore, hardness ratios of the folded lightcurve (see Fig. 6) do not 
show absorption (or hard X-ray excess) at the dips of the lightcurve, which implies 
that the flickering in X-rays is not associated with occultation of blobs. Our main 
conclusion is that the X-ray emission may be associated with an accreation disk 
corona rather then directly with the accreation column. In order to improve our 
understanding of the flickering oscillations, multiwavelength simultaneous obser- 
vations of TT Ari are required. 
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